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Enhanced Quantum Sensitivity in a Vibrating Diatomic Molecule due to Rotational 

Amendment 
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Quantum sensitivity is an important issue in the field of quantum metrology where sub-Planck 
scale structures play crucial role in the Heisenberg limited measurement. We investigate the meso¬ 
scopic superposition structures, particularly for well-known cat-like and compass-like states, in the 
rotating Morse system where sub-Planck scale structures originate in the dynamics of a suitably 
constructed SU(2) coherent state. A detail study of the sensitivity analysis reveals that rotational 
coupling in the vibrational wave packet can be used as a probe to enhance the sensitivity limit in 
a diatomic molecule. The maximum sensitivity limit is identified with the rotational amendment, 
and a quantitative measure of the angle of rotation for different rotational levels is also given. The 
correspondence of the numerical result with the angle of rotation is also delineated in phase-space 
Wigner representation. 
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I. INTRODUCTION 


Improvement in parameter estimation has often led to 
scientific breakthroughs and technological advancement. 
Recent advances in experimental techniques allow us ac¬ 
cess to unprecedented levels of control over quantum sys¬ 
tems. Quantum metrology is the field which exac^ 
deals with the fundamental limits to measurement [l|. 
To reach the ultimate sensitivity limit, one can repeat 
the measurement process N times and take the average 
over the outcomes. It reduces the error which scales as 
I/y/N, known as the standard quantum limit. This is 
the ultimate limit one can reach using classical proper¬ 
ties. To push this boundary, one needs the help of quan¬ 
tum properties. In quantum metrology, special states, 
such as entangled or squeezed states, have been employed 
for estimation of these parameters to beat the standard 
quantum limit Qi. In this case, the sensitivity can 
be enhanced \/N times and can reach the Heisenberg 
limit. On the other hand, the Planck scale executes a 
fundamental role in quantum mechanics. Phase-space 
quasiprobability distributions of certain quantum super¬ 
position states reveal structures on a scale that is smaller 
than the Planck dimension. The existence and impor¬ 
tance of these small structures (called sub-Planck struc¬ 
tures) were first pointed out in Ref. Q. These smallest 
interference (sub-Planck scale) structures play a crucial 
role in high-precision parameter estimation and Heisen¬ 
berg limited measurement. Recently, sub-Planck scale 
structures have drawn considerable attention and have 
been found in different situations [iQ. In all these 
studies, the sub-Planck interference phenomena appear 
with a suitable combination of appropriate superposi¬ 
tion of coherent states (CSs) [l^. In our earlier stud¬ 
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ies, we found the existence of sub-Planck structures in 
a molecular system and showed their decoherence sen¬ 
sitivity BO- It involved the vibrational motion of a 
diatomic molecule described by the evolution of a suit¬ 
able wave packet. Governed by a nonlinear energy spec¬ 
trum, the initial wave packet breaks into mesoscopic su¬ 
perpositions at a later time and gives rise to sub-Planck 
structures in phase space. 

In recent years, vibrational dynamics of diatomic 
molecules has gained importance due to its potential ap¬ 
plication in quantum computation. For example, ultra¬ 
fast Fourier transforms can be performed using a fem¬ 
tosecond laser-driven molecule [H, [13. High-precision 
molecular wave-packet interferometry has been used to 
read and write the amplitude and phase information of 
wave functions Eliil, which is a vital task for quantum 
information processing and the development of quantum 
gates. 

A key concern, however, is the effect of rotational cou¬ 
pling on the vibrational motion of diatomic molecules. 
Our goal is to study the sensitivity limit vis-d-vis sub- 
Planck structures in phase space by introducing the ro¬ 
tational coupling with the vibrational motion of a di¬ 
atomic molecule. Visualization of the rovibrational dy¬ 
namics needs a three-dimensional scenario [20l - l^ , hence 
it is difficult to observe its dynamics in phase space which 
will require a six-dimensional configuration. To unravel 
this intricacy, we recall an appropriate model, called 
the rotating Morse system, which can describe the rota¬ 
tional and vibrational coupling nicely in one-dimensional 
symmetry [2^ - 1^ . The energy eigenstates of a rotat¬ 
ing Morse system in phase space are elucidated in Ref. 
1^ . To the authors’ knowledge, there is no study about 
the wave-packet dynamics of this system in phase space. 
In this study, we choose a model [13 where the effec¬ 
tive potential becomes minimum around a certain equi¬ 
librium internuclear distance, which is a function of the 
rotational quantum number j. This can satisfactorily de¬ 
scribe the coupling between the two degrees of freedom, 
i.e., the rovibrational interplay. This coupling is also cap- 
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tured in mesoscopic superpositions states, such as cat-like 
states and compass-like states [13, . The most sensi¬ 

tive structures in phase space, called sub-Planck scale 
structures, are found to determine the sensitivity limit 
of a quantum state. We have chosen the example of 
an iodine molecule (h), which is a uniquely suited seed 
molecule for laser-induced fluorescence and an appropri¬ 
ate CS wave packet is constructed to see the system dy¬ 
namics. We show the amendment in vibrational wave- 
packet dynamics due to the presence of rotational cou¬ 
pling. Significant advances have been made in manipulat¬ 
ing and controlling rotational population in rovibrational 
wave packets by using shaped femtosecond pulses and 
wave-packet interference [3^ . To make the effect of rota¬ 
tional coupling transparent, we have considered a single 
rotational level for a rovibrational wave packet . The 
sub-Planck dimension in mesoscopic superposition struc¬ 
tures is found to vary with the rotational quantum num¬ 
ber j. Maximum sensitivity is achieved for a particular 
value of j and a scheme is provided to find the exact ori¬ 
entation of the corresponding system in phase space. Ad¬ 
ditionally, the corresponding phase-space Wigner distri¬ 
bution is numerically calculated and delineated in phase 
space to further verify the orientation and the structures 
in the maximum sensitive state. 

The paper is organized as follows. We present a brief 
overview of ID rotating Morse system and its validity. 
We construct the corresponding SU(2) CS wave packet 
to analyze the dynamics of the CS and to explain the 
effect of rovibrational coupling in configuration space. In 
Sec. Ill, we study the rotational sensitivity in a vibrat¬ 
ing diatomic molecule through mesoscopic superposition 
structures. Specifically, we have focussed on cat-like and 
compass-like states where the sensitivity issue is explored 
at the sub-Planck level. The scaling law is verified and 
the maximum sensitivity limit is achieved for rotational 
amendment. A quantitative measure of the angle of ro¬ 
tation in phase space is also depicted. Furthermore, a 
numerical study shows the phase Wigner distribution, 
which reveals a nice correspondence with the angle of ro¬ 
tation. Finally, we end up with some conclusions in Sec. 
IV. 


II. DYNAMICS OF COHERENT STATE IN 
ROTATING MORSE POTENTIAL 

We start with the effective potential 14// {r), known as 
the rotating Morse potential, 

Veffir) = (i) 

2/ir^ 

The first part describes the well known Morse potential, 
an appropriate model for vibrating diatomic molecule. 
D is the dissociation energy of the molecule, tq is the 
equilibrium internuclear separation, and /3 is the range 
parameter. The second part stands for the centrifugal 
contribution of rotation. A description of the system 


can be obtained with a modified equilibrium internuclear 
distance rj and a dissociation energy Dj [2^. Using a 
semianalytical method [l^, one can find 


U = ^0 



A 

P^rlD 


; Dj = D-A 



P^rlDj 


( 2 ) 


where A = . yvh define A,- = and expand 

the centrifugal term of Eq. m around r = rj. Keep¬ 
ing terms upto second order, the Schrodinger equation is 
solved to obtain the eigen functions of the rotating Morse 
system as 


V'nj(y) = (3) 

where the variable y: y = = 

(0 < y<oo) and u = n is the vibrational 

quantum number, L‘^{y) stands for associated Laguerre 

polynomial, s = ^{cq - Eyj)Xyc 2 , and Xj = 

Here the constants are expressed as cq = ^Ajb'j — SAjbj + 
Aj, Cl = {3Ajb^ — 2.Ajbj + uD)/u, C 2 = {3Ajb^ — Ajbj + 
u'^D)/u^, and bj = {j3rj)~^ are dependent on quantum 
number/. Defining A/ = A/, one obtains the constraint 


condition, 2s -|- 2n = 2A/ — I. Nnj is the normalization 


constant: Nnj = 


/3(2\j -2n-l)r{n+l) 


-, 1/2 


r (2 Aj —n) 

The rovibrational energy eigen values Enj turn out as 


A'n,/ — 2-—(n -|- 1/2) —T 1/2)^ + cq 


C2 


(4) 


It is worth pointing out that in the absence of rotation, 
Co = 0, Cl = C 2 = D, the system describes a vibrating 
diatomic molecule, z.e., the well known Morse potential. 

Alternatively, one can compute rj numerically by solv¬ 
ing the transcendental equation 


dVeffjr) 

dr 




(5) 


These two sets of r/’s are plotted in Fig. [IJa), which 
show very good agrement for / < 160 for I 2 molecule. It 
shows that the rotational motion increases the equilib¬ 
rium distance [Fig. UDa)] and decreases the dissociation 
energy [Fig. [Ub)]. Physically, in the presence of the ro¬ 
tational centrifugal force, the two constituent atoms of 
a diatomic molecule tend to settle at a larger distance 
and are more prone to dissociate, reducing the amount 
of energy required to make them independent. 

We construct a rovibrational wave packet of I 2 
molecule, which is a CS, dependent on particular ro¬ 
tational quantum number. Many theoretical and ex¬ 
perimental investigations have been carried out on this 
molecule, in particular, Zewail and co-workers investi¬ 
gated rovibrational wave packet dynamics in the well- 
characterized electronic H0+ state [^. Lohmiiller et al. 
[ 2 ^ discussed the pump-probe experiment of I 2 at room 















3 


(a) 


Dj 


0 50 100 150 200 250 300 0 50 100 150 200 250 300 

j j 

FIG. 1: (Color online) (a) The variation of Vj with the rotational quantum number j. It implies that the numerical values, 
obtained by solving the transcendental Eq. ([5]) (red filled circles) match nicely with the approximate values from Eq. ((2} (black 
circles). It starts to differ for higher values of j (> 160). The change of the dissociation energy of the effective potential with j. 
We have chosen I 2 molecule where /3 = 0.9605 a.u~^, reduced mass fi — 11.56 x 10"^ a.u., ro = 5.716 a.u., and D — 0.0198 a.u. 
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FIG. 2: (Color online) Effective potentials for j — 0 (solid 
line), j = 45 (dashed line) and j = 65 (dotted line) are de¬ 
picted. Inset: Zoom of the corresponding 10th energy lev¬ 
els, where the wave packets at t = 0.25 Trev are shown as 
dark filled (j = 0), dashed (j = 45) and light filled (j = 65) 
plots, respectively. The potential and the internuclear dis¬ 
tance are in atomic units. The corresponding parameter val¬ 
ues are /3 = 0.9605 a.u., reduced mass = 11.56 x lO'^ a.u., 
ro = 5.716 a.u., D = 0.0198 a.u. and a = 1.6. 


temperature and the detection of fractional revivals us¬ 
ing a full-dimensional quantum wave packet. Here, we 
consider an initial rovibrational wave packet which is 
centered around the 10th vibrational energy level with 
j = 45. Under the laser polarizations magic angle con¬ 
ditions [12, it takes into account the vibrational as well 
as the rotational motions. Once the bound states of the 
potential are included, the dynamical symmetry group 
becomes SU{2). For Morse system, the corresponding 


SU(2) generators are given (S^. In this case, we find 
^ d 


J+ = 

J_ = - 
Jo = 


, (2s-l)-f-s(2s-l)-A, 

[dy y 

—j—(2s “t“ 1) — —s(2s -f 1) -j- Xj 

dy y 


s — 1 


s + 1 


“ d y . 


( 6 ) 


Jo is the projection operator of the angular momentum 
m: m = n — Xj + 1/2. 

We obtain the SU(2) CS by operating the displacement 
operator exp(aJ+ —a* J_) on the highest bound state n', 
defined by J+ipn',j(y) = 0, where a is the CS parameter. 
Temporal evolution of the CS wave packet is given, in 
the eigen function basis, by 


^iy,t) = '^di^pn,jiy)e (7) 

n—0 


where the weighting coefficients are evaluated as 

. ^ (-«)»'-» r nMr(2A, -n) 1 ^ 

" (n'-n)! [n\T(2Xj-n')_ 


( 8 ) 


The presence of a nonlinear term in the energy expression 
leads to interesting phenomena, called fractional revivals, 
which occur at some specific instances between two full 
revivals [Hjllll- The short-time evolution displays a clas¬ 
sical periodicity. The classical and revival time periods 
are, respectively, given by 

27rA ■ 

Tci = X-^7-^, and Trev = 27rAj/c2. (9) 

ZCl — C2 /\j 


At fractional revival times {p/q)Trev (where p and q are 
mutually prime integers), the wave packet breaks into a 
number of subsidiary wave packets. For even (odd) val¬ 
ues of q, the wave packet breaks into q/2 (q) parts. In 
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the inset of Fig. [2 the 10th vibrational energy levels for 
different rotational numbers are zoomed and the rovibra- 
tional coupling effect is shown at t = 0.25 Tr-eu, when CS 
is split into two parts. For j = 0, the two parts are situ¬ 
ated at 5.3 and 6.48 a.u. (dark filled plot). For j = 45, 
they come close to each other, situated at 5.62 and 6.36 
a.u. respectively (dashed line). For a greater value of j 
(j =65), the position-space probability structure looks 
completely different and shows oscillatory structure (light 
filled plot). The interpretation lies in the fact that the 
two split CSs oscillate inside the potential well in a back- 
and-forth motion. In the first quarter of the oscillation, 
they approach each other, while in the next quarter, they 
recede. At halfway point of the oscillation, they are re¬ 
flected from the potential well with a phase change of tt 
and again become counterpropagating. For j = 65, they 
overlap each other in the course of their oscillation and 
produce the oscillatory ripples, clearly visible in the in¬ 
set of Fig. [5] A single interference ripple has dimension 
^ 0.1 a.u. or 5.3 picometers. Although, the experimen¬ 
tal observation of small quantum interference structures 
is very challenging, similar interference ripples in picome- 
ter scale were recently visualized experimentally for the 
I 2 molecule [H, 



Sensitivity by using classical action 


r. This Wgner function representation can reveal inter¬ 
esting mesoscopic superposition structures of the CS at 
different times. In addition to its positive regions,the 
Wigner function can also possess negative regions for 
nonclassical states. In the course of time evolution, one 
obtains the Schrodinger cat-like state at l/4th of the re¬ 
vival time. Four-way break up or the compass-like state 
emerges at l/8th of the revival time. Sub-Planck scale 
structures appear in the Wigner function at the inter¬ 
ference region of these mesoscopic superposition states. 
These structures are alternate tiles of maxima and min¬ 
ima. For symmetric potentials such as the harmonic os¬ 
cillator, these tiles are rectangular in shape and one can 
easily find the area of these structures by multiplying 
two side arms, by measuring the distances between the 
zeros of the Wigner function. However, for an asymmet¬ 
ric potential, it is quite nontrivial. We find the zeros of 
the Wigner function around a particular structure (either 
positive or negative) by projective plots of the Wigner 
function in both the conjugate coordinates. Then we per¬ 
form a set of measurements and, finally, take the average. 
The smallest sub-Planck structures are formed due to the 
superposition of off-diagonal superposition structures in 
a compass-like state. 

In addition to the above procedure, one can follow an 
alternative methodology, mentioned in one of our papers 
[HI- The idea is as follows: In typical experimental sit¬ 
uations, a small perturbation can be applied through a 
weak constant force, which will physically shift the state 
in phase space. This can be mathematically incorporated 
by finding an appropriate displacement operator for the 
coherent state, then applying the operator on the state 
for a small displacement. When the state is displaced 
by the length of a sub-Planck structure, the two states 
become quasiorthogonal and distinguishable. Hence, it 
decides the minimum amount of perturbation and force, 
which the system can detect. The overlap between the 
initial and final states in terms of the Wigner distribution 
is 


FIG. 3: Verification of scaling law between the sensitivity 
measured from Wigner plot and numerically calculated by 
evaluating the classical action. Here, the proportionality fac¬ 
tor is 3.78 and slope is 0.99('~ 1.0). 


III. ROTATIONAL SENSITIVITY 


Until now, we have explored the wave-packet dynam¬ 
ics in position space only. For a better description, we 
present a phase-space picture of the dynamics. Here, we 
make use of the Wigner function 37| , which is defined as 


1 r°° 

W{r,p,t) = ^ / <i)'*(r — r',t) 

7-00 

x<^'{r + r',t)e-^^P'^'/^dr'. ( 10 ) 


Here, <I>'(r, t) is the coherent state as a function of 


Q pOO POO 

($'(r,t)|$"(r,t)) = — / / W'{r,p,t)W''{r,p,t)drdp 

J — 00 J — 00 


( 11 ) 

where prime denotes the state before perturbation and 
double primes stand for the same after the application 
of an external perturbation. One can try to find the 
displaced state and its Wigner function after carrying out 
a lengthy calculation. The overlap function is oscillatory 
and the period of each oscillation will give the length of 
the structure in some particular direction in phase space. 
This direction is exploited by utilizing the complex form 
of the coherent-state parameters. One can, in principle, 
find the length in a number of phase-space directions to 
have a better idea of the shape of the structure. The 
method is cumbersome and not worthy to apply in the 
context of the present application because we need many 
of such kind of measurements. In the following, we will 
investigate how rotational coupling affects the quantum 
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FIG. 4: Variation of the sub-Planck dimension with j: (a) sub-Planck variation in cat-like state and (b) sub-Planck variation 
in compass-like state, (c) Comparison between these two cases. 



FIG. 5: (Color online) Points show the variation of all min¬ 
ima or most sensitive sub-Planck structures in cat-like and 
compass-like states. The solid line is drawn for cat-like state 
minima where alternate points depict majcimum sensitivity. 
Among them, second minima corresponding to j = 82 gives 
the most sensitive sub-Planck region in cat-like state. The 
dotted line shows the compass-like state minima variation 
where it brings out the most sensitivity for j = 64, occurs 
in the first minima. 


sensitivity. Specifically, we explore the sensitivity limit 
due to rotational amendment in the vibrating molecule. 


A. Mesoscopic superposition states and their 
sensitivity 

In this section, we make a quantitative analysis of 
the sensitivity of sub-Planck dimension with rotational 
coupling. Here, we denote the dimension of the small¬ 
est structure by ‘s’, which is propotional to ^ fi^/A 
1/A in atomic units), where A is the classical action 
of the state in phase space iini- The classical action 
is defined by the product of the effective support of its 
state in position and momentum spaces: A ^ Ax x Ap, 
where Ax = \J (x^) — and Ap = (jA) — (p)^. 
These quantities should be evaluated on the basis of the 
coherent-state wave packet given by Eq. 0. 

Now it brings out the question of evaluating the sen¬ 


sitivity or the area of the sub-Planck scale structures. 
In principle, the area can be estimated by either (i) 
measuring the area of the structures in the phase-space 
Wigner distribution or {ii) measuring the classical ac¬ 
tion. The first approach needs a huge computational 
time to evaluate the Wigner function integral in each 
case, upon choosing a proper phase-space region. Hence, 
it should be avoided, when one needs a large number of 
data. On the other hand, the latter approach requires 
the proper scaling law between the actual sub-Planck 
area and the classical action. In our technique, we have 
made use of both of the approaches to prevail over the 
situation. In the first step (i), we plot the Wigner func¬ 
tion at l/8th revival time for only six chosen j values 
(j = 0, 64,94,116,136,150) and measure the area of sub- 
Planck structure in each case as a reference value. In 
step (fi), we compute the inverse of classical actions for 
the same parameters and, finally, compare with the ref¬ 
erence values obtained from Wigner plot. The scaling is 
depicted in Fig. El which produces slope ~ 1.0. Hence, 
this is a confirmation of the scaling law: the inverse of the 
classical action is directly proportional to the quantum 
sensitivity. Now one would be able to perform a thor¬ 
ough quantitative estimate of the sensitivity. A system¬ 
atic analysis is performed to quantify the dependence of 
sensitivity with different rotational angular momentum 
quantum numbers and at different evolution times. 

Variation of sub-Planck dimension with the rotational 
coupling is shown in Fig. 0] As displayed in Fig. [1] 
here we keep increasing the value of j up to 160. Fig¬ 
ure EKa) shows the variation of the smallest interference 
region with rotational amendment for a cat-like state. 
Points depict the numerical values of sub-Planck dimen¬ 
sions which vary with the rotational coupling parameter 
j. It is interesting to see that the variation follows an 
oscillatory behavior where all the minima represent the 
high sensitive regions. With increasing value of j, all 
minima acquire comparatively higher values. In the cat¬ 
like state, the first minimum occurs at j = 38. It is 
noteworthy to mention that the presence of rotational 
coupling with particular j values corresponding to the 
minima shown in Fig. praise the sensitivity limit as com¬ 
pared to the case for j = 0. Figure EKb) gives the varia- 
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FIG. 6: (Color online) Time evolution of Wigner function of the CS. The first row shows cat-like state at t = for (a) 

j = 0, and (b) j — 82 and the second row shows compass-like state at time t — Trev/S for (c) j = 0, and (d) j = 64. The 
mesoscopic superposition of maximum sensitivity in Rotating Morse potential 


tion of sub-Planck region in a compass-like state. In this 
case, numerical data shows the oscillatory nature where 
all the minima capture the regions of greater sensitivity. 
The first minimum corresponds to j = 64. A comparison 
is made between these two cases in Fig. Uc). It shows 
that the compass-like state brings out the maximum sen¬ 
sitive state. A detailed study is given in Fig. [5l In the 
cat-like state, minima or the most sensitive sub-Planck 
dimensions are depicted by the points which are joined 
by a solid line. It shows that j = 82 brings out the most 
sensitive region in the cat-like state. In the compass-like 
state, minima are joined by a dotted line and it shows 
that the most sensitive sub-Planck region arises at the 
first minimum corresponding to j = 64. In the next sec¬ 
tion, further study involves the examination of the ori¬ 
entation of the system in phase space due to rotational 
coupling. 


TABLE I: Angle of rotation of the wave packet corresponding 
to the black dots in Fig. [5] for both cat-like and compass-like 
states. 


Cat State 

Compass State 

j 


j 

<t> 

38 

0.16 TT 

64 

0.22 TV 

82 

0.72 TT 

104 

1.16 TV 

116 

0.72 TV 

126 

1.71 TT 

142 

2.16 n 

150 

1.21 TV 

160 

2.77 TT 


B. Angle of rotation 

The rotational quantum number introduces rotation of 
the wave packet in phase space and in the above section 
we have found the states with maximum sensitivity for 
some particular values of the rotational quantum num¬ 
ber. Hence it is worth finding out the exact amount of 
rotation (/) corresponding to the states of maximum sen¬ 
sitivity. Here, we provide a numerical estimation of this 
rotation angle. It is well known that Jq is the generator 
of rotation and is related to the angular momentum [s^ . 
The corresponding rotation operator would be 17 = 

This operator upon operating on the initial wave packet 
gives 

n' 

n—0 

= x{y,i)- ( 12 ) 

The resulting state is rotated by an angle, depending 
implicitly on j. 

There is a one-to-one correspondence between the 
above state and the wave packet <!)(?/, t), directly obtained 
from the time evolution. Hence we find the angle of ro¬ 
tation by maximizing the overlap |(y(j/,t)|$(y,t))p for 
a given j. Numerically estimated angles of rotation for 
specific important values of j are shown in Table-H] The 
cat-like states reveal maximum sensitivity for j = 82 for 
which the rotation angle is found to be 0.72 tt. The an¬ 
gle for the most sensitive compass-like state (j = 64) is 
0.227r. 
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C. Phase space picture 

To obtain greater insight into what has been predicted 
in the previous section, we again invoke the phase-space 
picture. Figures [Bfa) and IHIb) display the wigner dis¬ 
tribution functions of the cat-like state for j = 0 and 
j = 82, respectively. Figure IHlJb) clearly shows rotation 
of the wave packet in phase space due to rovibrational 
coupling. It shows 0.727r rotation in anticlockwise direc¬ 
tion. Similarly, Figs. [6jc)-(d) show the wigner distribu¬ 
tion functions of the compass-like state for j = 0 and 
j = 64, respectively. Following the sensitivity study, we 
found that the compass-like state for j = 64 provides 
the maximum precision in this rotating Morse system. 
Although rotation of this particular state is obtained as 
0.227r from Table-Hl there is another crucial factor: the 
ratio of the revival and classical time scales. This ra¬ 
tio is not an exact integer in most cases and the extra 
fraction introduces an additional phase in the evolution. 
The extra rotation is calculated to be 0.0297r at 1/8 Trev 
When added to the rotation due to j, the resulting state 
is expected to rotate by 0.2497r ^ 1/4 tt, which is in 
conformity with the Wigner function in Fig. [S^d). 

IV. CONCLUSION 

Proper resource accounting is crucial when investigat¬ 
ing the precision or sensitivity in quantum systems and 


formulating the ultimate limits in quantum metrology. 
In this study, we have considered the rotational cou¬ 
pling in the vibrating diatomic molecule {I 2 ) and ex¬ 
plored the sensitivity of mesoscopic superposition struc¬ 
tures. Special attention is paid to cat-like and compass¬ 
like states where sub-Planck scale structures exist in the 
quantum interference region. Our sensitivity analysis of 
quantum interference structures reveals the fact that ro¬ 
tational coupling enhances the sensitivity limit in a vi¬ 
brating diatomic molecule. We have also identified the 
rotational levels corresponding to the maximum sensi¬ 
tivity limit. Our study avoids the complication of six¬ 
dimensional phase space for rovibrational dynamics of 
a diatomic molecule. The 1-D rotating Morse potential 
can well capture the rotational effect throughout the time 
evolution in phase space. Moreover, we provide a quan¬ 
titative measure of the angle of rotation for different an¬ 
gular momentum states. Our numerical result shows a 
nice correspondence between the angle of rotation and 
the phase-space Wigner representation. This study leads 
to an enhancement in the sensitivity limit and hence 
provides improvement in the Heisenberg limit for quan¬ 
tum metrology, which is not possible without rotational 
amendment. 
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